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ABSTRACT 


The  inductive  and  radiative  approaches  to  wireless  power  transmission  (WPT)  were 
simulated  using  commercial  software.  For  the  inductive  approach,  working  at  100  kHz, 
efficiencies  over  90%  were  obtained  at  short  ranges  (less  than  30  mm)  utilizing  ferrite  plates. 
For  the  radiative  approach,  the  transmission  loss  between  antennas  was  less  than  1  dB  at  ranges 
less  than  3  m  when  near  field  focusing  was  employed.  The  results  for  the  two  approaches  are 
important  because  they  demonstrate  that  efficient  transmission  of  energy  for  battery  charging  can 
take  place  between  a  WPT  ground  station  and  client. 
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I.  INTRODUCTION 


A.  BACKGROUND 

There  are  numerous  advantages  of  wireless  power  transfer  (WPT)  for  many  remote  energy 
source  and  battery  charging  applications.  The  approach  was  first  proposed  for  vehicle 
propulsion  in  mid  1960s  and  shortly  after  demonstrated  by  researchers  at  Raytheon  [1],  It  gained 
renewed  interest  in  the  early  1990s  for  micro-UAV  propulsion.  More  recently,  WPT  has  been 
used  for  charging  wireless  devices,  and  commercial  WPT  charging  products  have  appeared  on 
the  market  (for  example,  Witricity  and  Energous1). 

The  basic  block  diagram  of  a  WPT  system  is  shown  in  Figure  1.  Prime  power  is  provided  by 
the  base  station  (master),  converted  to  radio  frequency  (RF)  and  then  input  to  a  coil  or  antenna. 
On  the  client  side  the  received  power  is  conditioned  (filtered  and  transformed  in  voltage  and 
current)  and  subsequently  delivered  to  the  battery  or  power  plant. 


Base  Station 
(Master) 


Transmitting 


Rectenna 


Figure  1.  Block  diagram  of  a  WPT  system  for  battery  charging. 


Client 

(Slave) 


There  are  two  broad  approaches  to  the  problem  based  on  the  operating  range  (distance)  of 
interest.  They  are  the  inductive  and  far-field  (radiative)  regimes. 


B.  INDUCTIVE  APPROACH 

For  short  ranges  (less  than  a  couple  of  cm)  inductive  WPT  systems  are  viable.  Inductive 
systems  use  two  coils  (as  shown  in  Figure  2)  with  one  located  in  the  charging  station  and  the 
second  in  the  device.  Energy  is  transferred  by  the  magnetic  fields  linking  the  coils.  Essentially 
it  is  a  transformer  without  a  core  material.  At  the  receiving  coil,  circuits  are  required  to  rectify 
and  condition  the  output  voltage  for  charging  the  battery. 


See  http://www.witricitv.com/  and  http://energous.com/ 
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Inductive  systems  generally  operate  at  low  frequencies  (<  10  MHz  and  usually  in  the 
hundreds  of  kHz  range).  Efficiencies  greater  than  75%  are  routinely  achieved,  with  some 
systems  achieving  efficiencies  greater  than  95%,  however,  they  operate  at  very  short  ranges 
(several  mm).  The  short  range  can  limit  the  application  of  WPT  and  furthermore  alignment  and 
spacing  of  the  coils  is  critical. 


Voltage  induced  at 
terminals  of  client,  Vind 


Magnetic  flux  linked  to 
the  client  coil 


Number  of  turns  in  the 
master  coil 

Magnetic  flux  density  set  up 
by  the  master’s  current  I\ 

Figure  2:  Inductive  power  transfer  (after  [2]). 

The  transfer  of  power  is  governed  by  Faraday’s  Law  [2],  which  gives  the  induced  voltage  in 
the  client  coil  as 


V  = 

y  ind 


21 


dt 


(i) 


where  d>2l  is  the  magnetic  flux  set  up  by  the  current  in  the  master  coil  that  is  linked  to  the  client 
coil.  The  linked  magnetic  flux  is  a  complex  function  of  the  coil  geometries,  size,  frequency  of 
operation  f  and  distance.  Generally  it  is  not  possible  to  compute  021  in  closed  form  except  for 

some  very  simple  cases.  However,  electromagnetic  simulation  software  is  capable  of 
numerically  solving  for  the  currents  and  voltages  in  the  WPT  circuit. 


C.  RADIATIVE  APPROACH 


Radiative  WPT  systems  use  two  antennas  rather  than  coils,  and  the  energy  is  transferred  by  a 
propagating  wave,  as  depicted  in  Figure  3.  The  received  power  at  antenna  separation  d  is  given 
by  the  Friis  equation  [2] : 


_PtGtGrk\  .2 
r  (4  ndf  1  1 


(2) 


where  Pt  is  the  transmitted  power,  G,  the  charging  station  antenna  gain,  G,  the  receiving  (client) 
antenna  gain,  and  X=d  f  is  the  wavelength  ( c  is  the  speed  of  light, /is  the  frequency  in  Hz).  F  is 
a  factor  that  accounts  for  wave  propagation  effects  in  the  medium,  such  as  multipath  and 
attenuation.  L  is  a  system  loss  factor  that  includes  device  losses,  rectifier  efficiency,  etc.  The 
receive  side  the  antenna  has  an  integrated  rectifier,  and  is  called  a  rectenna. 

Radiative  systems  operate  at  higher  frequencies  than  inductive  systems  (>  1  GHz),  and  suffer 

a  1  Id1  propagation  loss.  Normally  d  must  be  great  enough  so  that  the  antennas  are  in  each 
other’s  far  field  region.  Higher  gain  antennas  can  be  used  to  increase  the  received  power,  but 
high  gain  antennas  have  a  farther  near-field  boundary,  and  they  become  physically  larger.  The 


2 


use  of  antennas  has  advantages  though.  Solid  state  array  antennas  allow  full  control  of  the 
antenna  excitation,  which  permits  scanning  and  focusing  of  the  beam.  This  capability  relaxes  the 
alignment  requirements  between  the  two  antennas,  and  can  be  employed  to  improve  the  system 
performance  when  the  antenna  far-field  criteria  are  not  met. 


Figure  3:  Example  of  a  radiative  WPT  system  (i.e.,  “microwave  beaming”  from  a  transmitting 
radar). 

A  disadvantage  of  radiative  systems  is  that  it  they  are  more  susceptible  to  losses  in  the 
medium  (i.e.  attenuation).  In  clear  air  there  is  no  significant  attenuation  for  frequencies  less  than 
15  GHz.  To  keep  the  losses  acceptable,  a  clear  line-of-sight  (LOS)  in  air  is  desirable.  Therefore, 
this  approach  cannot  be  used  for  vehicles  submerged  in  water  or  buried  in  wet  ground.  However 
is  can  be  used  for  ground  vehicles,  air  vehicles  on  the  ground,  or  even  warfighter  packed 
equipment.  It  does  not  require  the  precise  alignment  of  the  charging  station  and  client,  and  can 
be  designed  to  operate  at  distances  of  several  meters  or  more. 

D.  SAFETY  AND  INTERFERENCE 

Other  issues  that  must  be  considered  when  using  electromagnetic  energy  are  safety  and 
interference.  Radiative  systems  are  more  directive  than  inductive  systems  (significant  power  in 
the  main  beam;  negligible  power  outside  of  the  main  beam).  The  low  frequencies  used  by 
inductive  systems  penetrate  materials  more  readily,  including  the  human  body.  However, 
because  of  the  short  ranges  and  relatively  low  power  involved,  safety  should  not  be  an  issue. 
Both  types  of  systems  can  cause  interference  in  neighboring  electronic  systems.  Again,  because 
of  the  short  range  of  the  systems,  the  interference  introduced  by  a  practical  system  should  be 
limited  to  same  platform  (self)  interference. 
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II.  INDUCTIVE  WPT  STUDY 


A.  SYSTEM  DESIGN 

We  have  investigated  an  inductive  WPT  system  for  charging  underwater  vehicles  [3],  The 
concept  has  been  studied  previously  by  SSC  Pacific  in  San  Diego  [4],  It  is  a  short  range  (<  4 
cm)  system  that  operates  at  100  kHz.  The  low  frequency  is  dictated  by  the  requirement  that  the 
system  operates  efficiently  through  seawater. 

There  are  four  types  of  compensation  network  combinations  that  can  be  used:  1)  series- 
series,  2)  series-parallel,  3)  parallel-series,  and  4)  parallel-parallel.  Two  types  of  compensation 
networks  are  shown  in  Figure  4.  We  have  selected  the  series-series  approach  because  of  its 
simplicity,  and  it  is  the  topology  used  by  SSC  Pacific. 


Series-series 


Mutual 

Mutual  inductance 

inductance  (m) 

m 

Figure  4:  Equivalent  circuit  of  the  inductive  WPT  system  with  two  types  of  compensation 
networks. 


The  efficiency  is  the  power  delivered  to  the  load  relative  to  the  power  transmitted  [3]: 

pL  RMMf 

Pt  Rl(R2  +  Rl)2  +  (g>0M)2(R2  +  Rl) 


(3) 


where 

Vs ,  Voltage  supply, 

Rx ,  Internal  resistance  of  the  transmitting  coil, 

R2 ,  Internal  resistance  of  the  receiving  coil, 

Zj ,  Inductance  of  transmitting  coil, 

L2 ,  Inductance  of  the  receiving  coil, 

Cj ,  Compensating  capacitance  for  the  transmitting  coil, 

C2 ,  Compensating  capacitance  for  the  receiving  coil, 

/, ,  Transmitting  coil  current, 

I2 ,  Receiving  coil  current,  and 

M,  Mutual  inductance  between  the  transmitting  and  receiving  coils. 
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The  magnetic  coupling  coefficient  expresses  how  well  the  transmitting  and  receiving  coils 
are  magnetically  coupled  together.  The  higher  the  value  of  k  the  better  the  magnetic  coupling  is 
between  the  two  coils.  The  magnetic  coupling  coefficient  is  defined  as 


k 


M 

Vaa 


(4) 


We  can  show  that  there  is  an  ideal  matched  load  resistance  for  a  given  gap  (coil  spacing). 
The  resulting  RL  matched  can  be  expressed  as 


R  I., matched  ^ 


2+(aM)2^-. 


(5) 


From  Eq.  (5)  it  is  apparent  that  the  matched  load  changes  with  distance  (through  M).  Mutual 
coupling  M  (and  hence  k  which  depends  on  M)  are  not  easily  computed,  and  must  be  determined 
by  simulation  or  measurement. 

Several  commercial  software  packages  can  be  used  to  simulate  the  WPT  circuit  and  coils. 
CST  Microwave  Studio  was  selected.  The  frequency  domain  solver  was  used  to  find  the  voltage 
and  current  at  the  load  for  a  known  applied  voltage  at  the  input.  The  numerical  solver  provides  a 
rigorous  solution  of  Maxwell’s  equations  for  the  specified  boundary  conditions.  However  some 
numerical  errors  are  present  because  of  the  meshing  and  other  approximations  in  evaluating 
equations.  In  most  cases  these  errors  can  be  made  as  small  as  desired  at  the  cost  of  increased 
simulation  time  and  memory  requirements. 

The  MWS  model  is  shown  in  Figure  5(a).  The  coils  are  backed  by  plates  and  have  cores. 
The  plates  and  cores  are  shown  in  blue,  which  is  set  as  free  space  initially  to  model  the  coils  in 
air.  The  efficiency  is  given  by  [3] 


where  VL  =  voltage  across  the  load 
Vin  =  voltage  across  the  input 
Iin  =  current  at  the  input 

cos  0  =  power  factor  (cosine  of  phase  angle  of  voltage  minus  phase  angle  of  current). 

B.  SIMULATIONS  AND  MEASUREMENTS 

The  computed  efficiency  is  shown  in  Figure  5(b)  for  the  coils  in  air  with  a  spacing  of  16  mm 
(face  to  face).  The  computed  optimum  load  value  is  29  Q .  The  optimum  load  was  computed 
from  the  dimensions  of  the  coils  used  in  the  lab  measurements,  which  are  only  approximated  in 
the  MWS  model.  (It  is  not  possible  to  model  wires  in  MWS  as  thin  as  those  used  in  the  actual 
coils.)  From  the  simulation  data  the  optimum  appears  to  be  about  45  Q . 
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Lumped  element 
compensation  networks 


Simulated  Efficiency  vs  RL  at  16  mm 


Figure  5:  (a)  CST  MWS  model  of  the  WPT  coils  and  circuit,  and  (b)  simulated  efficiency  versus 
load  resistance. 


Measurements  were  performed  in  the  Microwave  Laboratory  by  Cena  [3]  and  the  results  are 
reported  in  his  thesis.  The  measurement  setup  is  shown  in  Figure  6(a)  and  the  results 
summarized  in  Figure  7.  At  16  mm  spacing  the  simulated  and  measured  data  are  in  close 
agreement,  but  both  are  slightly  below  the  calculated  value.  See  [3]  for  details  on  the  calculated 
value. 


Figure  6:  (a)  Measurement  setup  used  to  obtain  efficiency,  and  (b)  coil  installed  in  a  UAV  hull 
(from  [3]). 


Distance  between  the  transmitter  and  receiver  coils  (mm) 


-Measured 

-Calculated 


Figure  7:  Summary  of  efficiency  results  based  on  measured  data  (from  [3]). 
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A  way  to  confine  the  magnetic  flux  and  increase  coupling  is  to  add  a  ferrite  surface  behind 
the  receiving  coil.  The  placement  of  the  ferrite  board  behind  the  receiving  coil  was  suggested  by 
SSC  Pacific.  The  effect  of  a  ferrite  plate  is  illustrated  in  the  plots  shown  in  Figure  8.  The 
magnet,  which  is  equivalent  to  a  current  loop,  represents  the  transmitting  coil.  The  long 
rectangle  on  the  right  is  the  location  of  the  ferrite  tile  and  is  assigned  jur=\.  A  relative  magnetic 
permeability  of  1  represents  a  WPT  system  without  ferrite  tiles.  In  Figure  8(b),  the  long 
rectangle  on  the  right  is  assigned  //.  =  20  which  represents  a  system  with  ferrite  tiles.  The 

concentration  of  the  magnetic  field  in  the  red  area  (location  of  the  receiving  coil)  is  apparent 
when  comparing  the  two  figures.  The  plots  in  were  created  using  Vizimag  3.18. 


Place  receive  coil  in  With  magnetic 

this  area  of  material  in  this 

concentrated  flux  area 


No  magnetic 
material  in  this 


area 


Figure  8:  Application  of  ferrite  plates  to  concentrate  the  magnetic  flux  (after  [3]). 


MWS  simulations  were  performed  with  ferrite  plates  and  cores.  A  sample  of  the  results  is 
shown  in  Figure  9  as  a  function  of  load  resistance.  The  simulations  can  be  used  to  find  the 
optimum  load  resistance  from  the  peaks  of  the  efficiency  curves.  Note  that  the  efficiency  is 
improved  using  both  the  ferrite  plates  and  cores.  It  was  also  found  that  increasing  the 
permeability  of  the  tiles  improved  the  efficiency  and  reduced  the  efficiency’s  dependency  on  the 
load  value.  However,  a  small  efficiency  improvement  (e.g.,  less  than  10  %)  may  not  be  worth 
the  cost  in  terms  of  weight  and  complexity. 
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Ferrite  core 


(a)  (b) 

Figure  9:  Simulation  results  for  coils  with  ferrite  cores  and  plates. 

Measurements  were  taken  by  Cena  with  ferrite  plates  behind  the  receiving  antenna.  The  plate 
arrangement  and  a  summary  of  the  results  are  shown  in  Figure  10.  The  number  of  available  tiles 
was  limited,  which  forced  the  irregular  arrangement  shown  in  Figure  10(b).  The  measured  data 
also  showed  an  increase  in  efficiency  using  the  tiles  compared  to  air. 


Distance  between  the  coils  (mm) 


Calculated  power  transfer 
efficiency 


Measured  power  transfer 
efficiency  with  ferrites 


X  Measured  power  transfer 
efficiency  without  ferrites 


Measured  power  transfer 
efficiency  with  an  Al  UUV 
hull 


(a) 

Figure  10:  Measured  results  with  ferrite  cores  (from  [3]). 


(b) 


C.  SUMMARY  OF  INDUCTIVE  WPT 

The  measurements  and  simulations  of  inductive  WPT  have  shown  that  high  efficiencies  can 
be  obtained  at  close  coil  spacing.  A  frequency  of  100  kHz  was  used  to  allow  the  system  to 
operate  in  seawater  without  suffering  decreased  efficiency  due  to  the  water  resistance.  The  next 
phase  of  research  should  concentrate  on  designing  a  practical  interface  between  the  coils,  and 
optimization  of  the  rectifying  and  charging  circuit. 
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III.  RADIATIVE  WPT  STUDY 


A.  SYSTEM  DESIGN 

In  a  radiative  system,  the  master  station  transmits  a  propagating  electromagnetic  wave  to  the 
client.  Generally  the  master  station  has  a  high  gain  antenna  that  must  be  continuously  pointed  in 
the  direction  of  the  client.  This  type  of  WPT  has  been  proposed  for  numerous  applications  such 
as  space  based  power  (Figure  11),  battery  charging,  spacecraft  recharging  and  station  keeping, 
and  direct  propulsion  of  UAVs  and  hovering  airships.  The  client  antenna  is  usually  of  low  gain 
and  small  physical  size  because  it  operates  on  a  mobile  platform  (e.g.,  UAV). 


Figure  11:  Space  solar  power  (SSP)  concept  whereby  energy  collected  by  solar  panels  is 
transmitted  to  Earth  at  microwave  frequencies  (from  [5]). 

Equation  (2)  governs  the  received  power  versus  the  system  parameters.  The  antenna  gains 
are  given  by  [2] 

G  =  AjiAJA2  (7) 

where  A  =  Ac  ,  effective  area 

e  a 7 

A  =  aperture  area  (m2) 

ea  =  antenna  efficiency  (0  <efl<l) 

0n  ~  A/ L,  half-power  beamwidth  (HBPW)  in  radians  in  the  principal  plane  where  the 
dimension  is  L. 


B.  NEAR  FIELD  FORMULATION 
1.  Transmitting  Array 

At  close  distances  the  antennas  may  be  in  each  other’s  near  field.  In  the  near-field,  the  full 
gain  of  the  antennas  will  not  be  achieved  if  uniform  amplitude  and  phase  weighting  is  employed. 
If  the  master  station  antenna  is  an  array  with  full  amplitude  and  phase  control  at  the  element 
level,  then  near- field  focusing  can  be  used.  Figure  12  shows  the  near- field  operation.  For  a 
general  case,  let  the  array  focal  point  be  at  /y(xy,yy,zy)  and  the  observation  (field  point)  at 
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R(x,y,z ) .  (Note  that  here  R  is  not  a  resistance,  but  the  standard  spherical  polar  radial  variable.) 
The  array  is  located  in  the  x-z  plane  and  has  a  rectangular  grid  with  Nx  by  Nz  elements 
(1  <  m  <  Nx  and  1  <  n  <  Nz )  spaced  dx  by  dz .  The  normal  vector  to  the  array  is  y . 


Figure  12.  Planar  array  antenna  near  field  model. 

Other  parameters  of  interest  are: 

•  The  position  vector  to  element  m,n :  dmn  =  xmnx  +  zmnz  where  xmn  = 

_2n-(Nz+l)A  _A  ^  _A 

zmn  2  ’  an<^  y nm  ^  • 

•  Position  vector  to  the  focal  point:  iy  =  xyx  +  yy  v  +  zj-z . 

•  Position  vector  to  the  observation  point:  R  =  xx  +  yy  +  zz . 

•  Vector  from  element  m,n  to  the  focal  point: 

R mn  ~  (xf  ~  xmn  (f  /  —  izf  ~  zmn  ) z  • 

•  Vector  from  element  m,n  to  the  field  point: 

=  (x-  xmn  )x  +  (y-  ymn  )y  +  (z-  zmn  )z 

so  that  tmn=5*n- 


2m  -  (Nx  + 1) 


d 


X  ’ 


(8) 


(9) 


Let  the  elements  be  z-directed  and  linearly  polarized.  A  general  form  for  the  element  factor 


is 


Fmn  (@mn »  fit 'mn )  z  ^ o  ’  ^  (rmn  )  ’  fmn  (@mn » 0mn  ) 


(10) 


where: 


fmn  ( 0mn  ,  fnn )  is  the  normalized  element  pattern  as  a  function  of  angles  6mn ,  <f>nm  . 
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q  j Prmn 

•  x¥(rmn)  = -  is  the  free  space  Green’s  function  with  J3  -  2n  /  A, . 

rmn 

•  C0  a  complex  constant  that  depends  on  the  element  type  (e.g.,  dipole,  patch,  etc.). 

•  0mn  is  the  angle  between  the  element  m,n  axis  (along  z)  and  the  field  point  position 
vector: 

COS  0mn  —  Z  D  fmn ,  sin  0mn  —  -^1  —  cos  (11) 


(j)mn  is  the  element  azimuth  angle: 

tan  (j)mn 


y f  ymn 

xf  —  xmn 


(12) 


To  this  point  we  have  allowed  for  the  patterns  of  each  element  to  be  different.  Typically  all 
elements  are  the  same  so  that  Fmn  ( 6mn ,<pmn)  =  F  {dmn, (f>mn )  or  fmn  ( 6mn ,<t>mn)  =  f(0mn, <Pmn )  • 

If  the  elements  are  unit  current  amplitude  z-directed  half-wave  dipoles  at  height  h  above  a 
ground  plane  that  is  a  perfect  electrical  conductor  (PEC),  then  the  ground  plane  element  factor  is 


7^0  z- 


J Prm 


COS 


71  a 
—  cos  0, 


mn 


v 

rmn 


sin  0, 


-  sinf  /?/?  sin  6mn  sin4„). 


(13) 


mn 


KVUmn) 


f(Fmn  '■tfinui ) 


In  (13)  the  far-field  condition  has  been  assumed  for  the  ground  plane  factor.  This  is  an  idealized 
result  for  an  isolated  element  above  an  infinite  ground  plane.  The  ground  plane  factor  for  a 
dipole  over  a  finite  ground  plane,  or  in  an  array  environment,  or  in  the  near  field  will  differ  from 
sin(/?/i  sin  6mn  sin  <f>mn )  due  to  mutual  coupling  and  diffraction.  These  effects  can  be  accounted 
for  by  modifying  the  ground  plane  factor  as  described  in  the  Appendix. 

The  total  array  field  at  the  observation  point  is  given  by  the  weighted  sum 

Nx  Nz 

E(R,e,t)= 

m= 1 n= 1 

'j  (14) 


NXNZ 

yi20fz  i- 

m= 1  n= 1 


o~jPr, 


where 

•  Amn  is  the  element’s  current  amplitude  coefficient. 

•  ®mn  Is  l'ie  phase  required  to  focus  the  beam  at  /y : 

^ mn  ~  f  ~rmn)  •  (15) 


•  y/mn  is  a  miscellaneous  calibration  and  correction  phase. 

In  Figures  13  and  14  we  plot  the  electric  field  intensity  for  an  array  of  half-wave  dipoles 
above  a  ground  plane  with  the  parameters  in  Table  1. 
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Table  1.  Parameters  for  the  plots  in  Figures  13  and  14  (extracted  from  Matlab  script) 

f=10e9;  %  frequency 

wave=3e8/ f ; 

%  Observation  parameters 
Xmax=2;  %  -Xmax  to  Xmax 

Ymax=6;  %  DZ  to  Zmax  to  avoid  singularity 

zp=0;  %  observation  plane  height 

%  Array  Parameters 

Nx=30 ;  Nz=2  0 ; 

dx=0.5*wave;  dz=0.5*wave; 

hgp=0 . 25*wave ; 


The  subtle  variations  in  the  field  distribution  are  more  easily  seen  if  the  free  space 
attenuation  of  the  field  is  removed.  Thus  at  an  observation  point  a  distance  R  from  the  origin 

gnonnY.pttt  MindB:  201°8l0  ^  | .  (16) 

In  Figure  13  is  shown  a  plot  of  the  normalized  electric  field  intensity  in  V/m  and  dBV/m  when 
the  transmit  array  is  focused  at  2  m.  In  Figure  14  the  result  shown  is  for  the  focal  point  at  2000 
m,  which  is  in  the  array’s  far  field  (i.e.,  a  linear  phase  is  used  to  scan  the  beam).  The  black  stem 
is  located  at  the  focal  point.  It  is  clear  that  the  field  reaches  a  maximum  at  the  focal  point  in 
Figure  13,  whereas  that  is  not  the  case  in  Figure  14.  There  is  a  3.3  dB  increase  in  signal  at  the 
focal  point  compared  to  far-field  scanning. 


(a)  (b) 

Figure  13.  (a)  Plots  of  normalized  electric  field  in  V/m  and  (b)  in  dBV/m  when  the  array  is 
focused  at  2  m. 
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(a)  (b) 

Figure  14.  (a)  Plot  of  normalized  electric  field  in  V/m  and  (b)  in  dBV/m  when  the  array  is 
focused  at  z  =2000  m  (i.e.,  the  far  field).  The  MR  factor  has  been  removed. 

2.  Receiving  Array 


Let  a  receiving  array  be  centered  on  the  y  axis  at  yp  with  numbers  of  elements  in  the  x  and 
z  directions  of  N'x  by  N'z  and  spacings  d'x  and  d'z  .  Similar  to  the  transmitting  array  the 


2p-(N'x+l) 


2q-(N'z  +  l) 


element  locations  are  given  by  x'pq=  — — v  x  d'x  ,  z'pq  =  — — z  —d'z  and  y'r)a  =  0 


rpq 


(l<p<N'x,l<q<N'z  )  as  shown  in  Figure  15.  The  primed  quantities  and  indices  p  and  q  are 
associated  with  the  receive  array. 


z  1 


Figure  15.  Total  field  at  element  p,q  of  a  receiving  array. 
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The  incident  field  at  element  p,q  is  obtained  by  setting  its  location  as  the  field  point  in  Figure 


12: 


E  (R,0,0)  — »  E  (x'pq,y'pq,z'pq)  =  Epq  . 


The  total  field  at  element  p,q  due  to  all  elements  in  the  transmitting  array  is  given  by  Eq.  (14) 
with  the  appropriate  angles  and  distances  inserted 


Nx»y 


Epq  =  yl20z  I  I 

m= 1 n= 1 


„  J Prmnpq  .  ¥ 

6  —  '  J^mnJWn 


'  J i^mnpq  ’  fimnpq  )  e 


rmnpq 


where 


( xpq  xmn)^  +  (y pq  ymn)9  +  (z pq  zmn)z  • 


'  mnpq  \Jvpq 


pq 


(17) 


(18) 


'mnpq 


mnpq 


and  0mnpq  is  the  angle  of  the  axis  of  transmitting  element  m,n  with  the  LOS 


^mnpq  ^z 


vector  to  receive  element  p,q  (an  extension  of  Eq.  (11)) 

COS  [^mnpq  )  — 

Similarly  extending  Eq.  (12),  the  azimuth  angle  is 

^„Mypq~ym" 


mnpq 


(19) 


xpq  xmn 


(20) 


The  complex  voltage  out  of  receiving  element  p,q  is  given  by 

Epq  E pq  D  hpq  ^  Omnpq ,  (f>mnpq  j  (21) 

where  hpq  is  the  vector  effective  height  of  receive  element  p,q  and  ( O'mnpq  >  fimnpq )  the  direction 
from  the  receive  element  axis  with  the  LOS  to  transmit  element  m,n 


cos(d'mnpq)  = 


r  Uz 
'mnpq^^  pq 


mnpq 


,  /  n  \  ymn  y pq 

tgn  (Knpq)  =  -,  • 


(22) 


(23) 


xmn  xpq 


Note  that  the  axes  of  the  receiving  array  (primed)  are  related  to  those  of  the  transmitting  array  by 

z  =  z,  x  =  —x,  y  =  —y  (24) 

and  for  every  point 

z'  =  z,x'  =  -V,  y'  =  -y .  (25) 
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For  a  half-wave  dipole  on  the  z'  axis  over  a  PEC  ground,  the  effective  height  is  derived  in 
the  Appendix: 


hpq  J  p  f  (  @mnpq  ’  fimnpq  )  @  • 


(26) 


This  is  an  idealized  result  for  an  isolated  element  above  an  infinite  ground  plane.  The  effective 
height  of  a  dipole  over  a  finite  ground  plane  or  in  an  array  environment  will  differ  from  this  due 
to  mutual  coupling  and  diffraction. 

The  final  step  is  to  calculate  the  total  complex  voltage  at  the  output  of  the  receive  array, 
which  is  the  complex  sum  of  all  receiving  dipole  voltages 


K  K 

Kut  =  Z  Z  V pq  A pq? 

p=\  q= 1 


,J¥pq 


480  Ay  ^A  Ay  Nz  (  /v  /  \ 

=  ;  — I  Z  Z  i:\zUffpq  hpq  y^mnpq  >  fimnpq  )  ^ pqe 


P  p= 1  q= 1  m=  1  «=1 


J¥ 


pq 


(27) 


x  ^( rmnpq  ')f^mnva  ’  ^mnva 


mnpq  >  rmnpq  '  inn 


J®mn  , mn  \ 


where  A'pq  and  y/'pq  are  miscellaneous  calibration  and  correction  amplitudes  and  phases  applied 

on  receive,  if  any.  0'pq  is  the  unit  vector  in  the  spherical  system  centered  at  p,q.  Given  that  Eqs. 
(23)  hold 

2  ^'pq  =  —  sin  O'mnpq  ■  (28) 

The  final  result  for  both  arrays  with  half-wave  dipoles  above  a  PEC  ground  plane  is 


480  K  K  „  jYpq 


^out  =  j  —pr  Z  Z  A'  e 

P  p=lq=\ 


NXNZ 


Z  Z  sin  O' 


m= 1  n=\ 


mnpq 


—  /  Br 

g  Jr' mnpq 

c°s(^wmnpq) 

rmnpq 

V 

/Tv 

1  vvmnpq 

sin  (Phvmnpq) 


v  a  J®  mn  JWmn 
^ mn e  e 


COS  I  ^Wm„pq 


l-(w') 


/\2 


sin  (phv'mnpq) 


mnpq 


(29) 


where 


mnpq 


=  sin  6 \ 


mnpq 


sin^ 


mnpq 


Wi 


mnpq 


—  COS  Omnpq  ?  sin  0mnr>a  =\ll~w> 


mnpq 


vmnpq 


mnpq 


sin  9mnpq  sin  (f>n 


rmnpq 


mnpq 


=  COS  O'mnpq  >  sin0'  =A H™') 


> mnpq 


(30) 
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Note  that  in  Eq.  (29)  the  voltages  of  the  elements  have  been  summed  (i.e.,  series 
combination).  For  parallel  combination  the  voltage  is  constant  and  we  should  use  ^out  KN'xK). 
The  time-average  power  is 


Pl=-Rc{VI 


2(N'sN'.) 


Re 


lit lit 


(31) 


The  power  delivered  to  a  conjugate  matched  receiving  load  of  resistance  RL  is  [6] 

p  _J!mL 

1  (KK)Ml' 


(32) 


C.  SIMULATION  RESULTS 


Several  benchmark  cases  were  run  to  validate  the  software  written  in  Matlab.  For  the 
validation  purposes  only,  a  frequency  of  300  MHz  is  used.  The  commercial  software  package 
Savant  was  used  to  compute  the  power  transferred  between  a  10  by  10  transmitting  array  and  the 
receiving  antenna  (see  Figure  16). 


Figure  16.  Savant  model  of  the  transmitting  array  and  receiving  antenna. 

A  comparison  of  the  results  is  shown  for  ranges  greater  than  10  m  in  Figure  17(a).  Equation 
(2),  the  Matlab  results  for  Eqs.  (29)  -  (32),  and  the  Savant  simulation  are  all  in  good  agreement. 
The  far  field  condition  is  assumed  to  hold  -  there  are  no  phase  corrections  for  focusing.  A 
comparison  of  the  focused  and  unfocused  cases  is  shown  in  Figure  17(b).  The  focused  case 
gives  an  improvement  of  12  dB  or  more  at  ranges  less  than  3  m.  This  confirms  the  results 
presented  previously  in  Figure  13. 
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(b) 

Figure  17.  Power  received  by  a  single  dipole  over  a  ground  plane  from  a  10  by  10  array  with 
ground  plane  transmitting  1  W.  (a)  Comparison  of  summation  formula  to  Friis  and  Savant 
simulations  (no  focusing),  (b)  Comparison  of  summation  formula  to  Savant  simulations  for  near 
field  ranges  with  and  without  focusing. 
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IV.  SUMMARY,  CONCLUSIONS  AND  FUTURE  WORK 


The  inductive  and  radiative  approaches  to  WPT  were  simulated  using  commercial  software. 
For  the  inductive  approach,  working  at  100  kHz,  efficiencies  over  90%  were  achieved  at  short 
ranges  (less  than  30  mm)  utilizing  ferrite  plates.  For  the  radiative  approach,  the  transmission 
loss  between  antennas  was  less  than  1  dB  at  ranges  less  than  3  m  when  near  field  focusing  was 
employed.  The  results  for  the  two  approaches  are  important  because  they  demonstrate  that 
efficient  transmission  of  energy  can  take  place  between  the  WPT  ground  station  and  client. 

The  frequency  used  in  the  simulations  was  300  MHz,  however  the  efficiency  conclusions 
apply  to  other  frequencies  when  the  antennas  are  scaled  in  wavelengths.  For  example,  the  10  by 
10  array  with  ground  plane  is  approximately  5  m  square  at  300  MHz,  but  would  reduce  to  0.5  m 
square  at  3  GHz,  or  0.15  m  square  at  10  GHz. 

The  next  step  in  the  research  (Phase  II)  is  to  design  efficient  rectifying  and  battery  charging 
circuits.  Candidate  designs  can  be  simulated  using  the  commercial  software  package  Advanced 
Design  System  (ADS)  from  Keysight  (formerly  Agilent).  Rectifier  efficiencies  of  90%  have 
been  reported  at  low  frequencies  (<  10  MHz)  because  efficient  diodes  are  available  in  this  range. 
However,  at  higher  frequencies,  (e.g.,  10  GHz),  achieving  this  level  of  efficiency  is  more 
difficult  [7,  8]. 
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APPENDIX  A.  EFFECTIVE  HEIGHT  OF  A  HALF-WAVE  DIPOLE  OVER 

A  GROUND  PLANE 

Consider  element  m,n  to  be  a  z  directed  dipole  at  a  height  h  above  a  ground  plane  that  is 
located  in  the  x-z  plane.  The  ground  plane  factor  is 


GPF  =  e 


jf5h  cos  a 


ymn  _ 


-jfih  cos  or 


ymn 


=  2/sin(/?/zcosai;  ) 

J  ymn  ' 


where  cos  ay  =  vmn  =  sin  6mn  sin  <j)mn  is  the  y  direction  cosine.  The  element  factor  for  a  z 
dipole  over  the  ground  plane  is 


EFmn  =  f2/ sin(/^ sin 6mn  sin <t>mn ) 


g  j Prmn 

=  -120  z - sin (Phvmn) 

rmn 


j6  0- 


o  j Prm 


COS 


f  ji  ^  ^ 

—  cos  ft 
2  7 


j 


cos 


'mn 


f  7t  ^ 

wmi 

V 2  j 


sin#, 


mn 


v 


1  -w' 


mn 


The  effective  height  is  determined  from 


v  =  f  uh 

r  mn  J^mnLJnmn 


and 


E  =  -  in  I—  h 

^ mn  J  /  rimn  • 


47rr 


lmn  ■ 


For  a  half  wave  dipole  with  a  ground  plane: 


cos 


h  —  —  / 

nmn  J 


r  n  ^ 

— cos  0mn 

V2  j 


P  sin  01 


sin(/3h sin Gmn  sin </>mn)0mn. 


mn 


Typically  we  choose  h  =  A  /  4  so  that  fth  =  n  1 2  .  This  is  an  idealized  far  field  result  for  an 
isolated  element  above  an  infinite  ground  plane.  The  effective  height  of  a  dipole  over  a  finite 
ground  plane,  in  an  array  environment,  or  in  the  near  field  is  different  from 
sin(/?/?  sin  6mn  sin  <f>mn )  due  to  mutual  coupling  and  diffraction.  These  effects  can  be  accounted 

for  by  modifying  the  ground  plane  factor  to  include  an  efficiency  factor  (A)  and  taper  exponent 
(. B )  to  match  the  actual  behavior: 


D 

A  sin  (fih  sin  6mn  sin  (j)mn ) . 
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